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^ , It is shown that the fuU unknown state of a spin-i system, S, which, within Horn's 

CO i statistical interpretation, is meant as the state of an ensemble of identically prepared 

' systems and is described by its density matrix, can be determined with a simultaneous 

' measurement with the help of an "assistant" system A whose initial state is known. 

, The idea is to let S and A interact with each other in a known way during a proper 

C " 3 , interaction time r, to measure simultaneously two observables, one of S and one of 

' A and their correlation. One thus determines the three unknown components of the 

^ ~ ' polarization vector of S by means of repeated experiments using a unique setting. In this 

""^5 1 wS'J can measure simultaneously all the non-commutative observables of S, which 

■ might seem prohibited in quantum mechanics. 

M-i. 

.^.^ ' Keywords: State determination; Quantum measurement; Two-level system; Coherent 

^ ' state. 

^3 : 

CT^' 1- Introduction 

^ ^ ' The determination of the unknown state of a quantum system is one of the most im- 

. portant issues in the field of quantum information.^'^ This determination involves 

\ a measurement process in which a macroscopic system, apparatus, is coupled to 

- - ' the quantum system; during this process the state of both the apparatus and the 

system is modified.'^ For instance, as currently described in many textbooks, the 
2-component of the polarization vector of a spin-^system, S, can statistically be 
determined by means of a repeated Stern-Gerlach experiment. In this process, the 
X- and ^/-component of the polarization vector are destroyed as a consequence of 
the non-commutation of the spin operator in the transversal directions. Other ex- 
perimental settings seem therefore necessary to measure the unknown polarization 
vector of S. Its three components are represented by incompatible observables, the 
Pauli operators, and their direct determination requires three macroscopic appara- 
tuses, differing by a change of orientation of the magnets and detectors. Likewise, 
the state of any two-level system, represented by a 2 x 2 density matrix p can be 
fully determined only through measurement of three linearly independent observ- 
ables which do not commute and cannot be simultaneously measured. 
Nevertheless, we will prove that the whole unknown density matrix of such a system 
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S, in particular the full polarization vector of a spin-^ system, can be determined 
indirectly by means of a set of measurements performed simultaneously on S and 
an auxiliary system, A, which we term the assistant. The strategy is the following: 
initially S is in the unknown state that we wish to determine, while the assistant 
A is in some known state. During some time interval S and A interact in a known 
fashion. Their joint state is modified, involves correlations and keeps memory of the 
initial state of S. A simultaneous measurement of one of the observables of S and A 
is then performed. Repeating this process provides then three statistical data: the 
expectation values of these observables and their correlation. We will show that one 
can infer the three components of the initial polarization vector of S from the three 
data. 

There are two approaches to this problem. Either using another two- level assistant, 
for solid state applications^ or using an electromagnetic field as an assistant, for 
quantum optics applications. 

After defining the general problem, we illustrate, in section 3, the first approach, 
i.e. using another spin-^ system as an assistant in a known pure state to determine 
the initial state of S. Then, in section 4, we show that it is also possible to use 
an electromagnetic field in a coherent state to determine the whole elements of the 
unknown density matrix of S. 



2. Statement of the Problem 

The idea of mapping the state p of an unknown spin-i system, S, onto a single 
observable of S+A system by using another system in a known state R, A, was first 
proposed by D'Ariano.'* It was explicitly implemented in a dynamic form in Ref. 5. 
The state of the composite system S + A, which is tested, is 

= Un^U^ (1) 

where the initial state of S* + A is TZa = R ® p and the evolution operator is 
U = e~'^^'^. Therefore, the dynamics of the system yields the required mixing of p 
and R and the simplest possible non-degenerate observable of the composite system 
S-l-A, O, can be chosen as a factorized quantity 

f) = w 6, (2) 

where the observable a) and 6 pertain to S and A respectively. Then the spectral 
decomposition of w and 6 can be used to construct the projection operator of O 

m n 
O) = ^ OJiTti, O = ^ OaPa, (3) 

i=l a=l 

where Pa and ni are eigen projectors of the observables 6 and u) respectively. There- 
fore, projection operator with a = {i, a) takes the form 

Pa = Pia = T^i ^ Pa. (4) 
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Repeated measurements of f2 which means repeated simultaneous measurements of 
Cj and 6, determines the joint probabihties to observe LOi for S and Oa for A 

F„ = P„ = Tr[7t,(7r,(K)pa)], (5) 

where TZr is defined in (Eq. [T]) . In fact the numbers Pa are the diagonal elements 
of tj\p R)U in the factorizcd basis which diagonalizes uj and o. 
The whole elements of the density matrix p can be determined by the mapping 
p ^ Pa- Like the idea of finding a universal observable,'' if H couples S and A 
properly, this mapping will be expected to be invertible for n > m. We shall see 
that even simple interactions can achieve this condition. For given observables cD of 
5* and d oi A and for a known initial state R, the precision of this procedure relies on 
the ratio between the experimental uncertainty of Pa and the resulting uncertainty 
on p, which can be characterized by the determinant. A, of the transformation 
([5]). For A = it is impossible to determine p from Pa- This means, the system is 
unstable with respect to small errors made during experimental determination of Pa 
or equivalently ((<t^) , (sz) , (szO'z))- Therefore, the Hamiltonian H and time interval 
T should be chosen so as to maximize |A| over all possible unitary transformations. 

3. Spin-^ Assistant in a Known Pure State 

In this section we illustrate the above ideas by studying a two-level system S, 
namely, an spin-i system in a known pure state. The density matrix of a spin-^ can 
be represented with the help of Pauli matrices. The determination of p corresponds 
to determination of the elements of the polarization vector, p- we let S and A 
interact during the time interval r. The observables a) and 6 to be measured are 
the 2-components of spin of S and A and are determined by az and -Sz respectively. 
The projection operators are 

■^i = ^(i + ^z), Pa = ^{i+Sz), (6) 

for i and a equal to ±1. Experiments will determine the four joint probabilities Pa = 

P++, P^ , P I-, P These probabilities are related to the three real parameters p 

of p by inserting (Eq. [6]) and R= i(l + s^) into (Eq. [5]). 

Pa = Ua + Va ■ P, (7) 

where 

Ua^^[U{i^R)U%.a, Va = ^[U0^R)U%,a, (8) 

with a — {ia} = {H — h, H — , — h, } and matrix elements have been represented 

in the standard representation of the Pauli matrices, & and s. The probabilities Pa 
should be positive and normalized for any the density matrix, p, such that < 1. 
These conditions imply that 

Ua>\Va\, ^Uq = 1, ^Va^O. (9) 
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The determinant of transformation p can be either 

«++ • X (10) 

or any other permutations of three of the vectors v^j^ , , v y. and v Therefore, 

the determinant of the transformation is four times the volume of the paraUelepiped 
made by three of these vectors. For example, 

A = 4i;++.(t;+_ xt;+_). (ii) 

If the unitary evolution operator U is such that vectors Va are not coplanar, the 
transformation (Eq. [7|) is invertible and one can determine p from the set of ■ 
Alternatively, p is deduced from {oz), (sz) and {sz(Jz) at time t. The notation SzC^z 
is used for simplicity instead oi Sz ® Oz which lives in the common Hilbert space 
of S and A. cTz, Sz and SzCfz can be simultaneously measured and are in one to one 
correspondence with the set of probabilities Pa ■ 

We first look for the upper bound of the determinant of transformation (Eq. [T]), 
|A| implied by the conditions (Eq. [9]). First we note that |A| increases with \va\ for 
each a. We therefore maximize under the constraints 

= ^z;„ = o. (12) 

This yields a symmetric solution for all these vectors 

Ua^\Va\ = -, COS ^ = --. (13) 

4 li'ttlh/jl 3 

By definition this means that vectors Va form a regular tetrahedron. These solutions 
are not unique and they follow from one another by rotating in the space of the 
spins and permutations of the indexes a. Therefore, the corresponding determinant 
for the upper bound is 

Having a non-zero determinant for the proposed procedure, ensures its feasibility. 
One simple choice for vectors Va is 

"-^ its"' '-^i^'-''-'''*' 

This yields a simple form for the density matrix of S: 

pi = \/3((j^), p2^Vi{sz), Pi = Vi{szcrz), (16) 

which gives directly the whole elements of the density matrix, /), in terms of the 
expectation values and the correlation of the commuting observables and Sz in 
the final state. 

Next step is to find out the interaction Hamiltonian and the interaction time r 
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which give such a description of the tested system, S. 

This correspondence can be achieved under the action of the Hamiltonian 

H = -^a-j:{sx cos(f> + Sz sincf) + ^[(sj^ - s^) sinc^ + coscj)], (17) 

where 20 is the angle between and the 2-axis, that is, cos(/) — Noting 

that = SIT? Xi where x satisfies cosx — 1/2 cost/), and taking as duration of the 
evolution r = x/sinx, we obtain U = exp^—iHr) — cosx — iH. The simpler form 

H = —^a^s^ + -(SySin(j) + s^) (18) 

of H can be obtained by a rotation of s and also achieves an optimal mapping 
p ^ Pa, provided Sz — > Sx smcf) + coscj) both in the measured projections Pa = 
1/2(1 ± s^) and in the initial state p+. The first term in (Eq. [T8|) describes, in 
the spin language, an Ising coupling, while the second term represents a transverse 
magnetic field acting on the assistant A. 



4. Assistant System as a Coherent State of Light 

In this section we discuss the possibility of using light as an assistant to determine 
the elements of the density matrix of a spin-i system. We show that in case of using 
an electromagnetic field in coherent state, one can determine the state of S from 
the commutative measurements on S and A. 

To describe this physical situation wc choose the Jayncs-Cummings model^'^ a 
widely accepted model describing the interaction of matter (two-level atom or spin- 
i system) and a single mode of radiation. This model is exactly solvable but still 
rather non-trivial and it finds direct experimental realization in quantum optics. 
The Hamiltonian reads :^ 

H = Ha + Hs + HsA = hoja^a + ^huaz + h'j{a+a + (T_a^) (19) 

where and a are the standard photon creation and annihilation operators of the 
field (the assistant A), with commutation relation [a, a^] = 1, ct^ are the standard 
Pauli matrices for the spin of the two level system S. The total Hamiltonian (|19p 
is the sum of the Hamiltonian of the field Ha, the Hamiltonian Hs of the two level 
system 5* and the interaction Hamiltonian Hsa which can be written as h'jV, with 
7 the coupling constant. It can be easily checked that the interaction operator, V, 
and the total number of excitations, N: 

V = d-+a + a-^a^ , N = a)a + a+a^. (20) 

are two constants of motion. V and N also commute with each other since = N . 
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One can then calculate exactly the relevant observables at time t using the 
Heisenberg equations of motion: 



a = —loja — ij(7-, 

cFz = 2i7(a^(T_ — (T+d). 



(21) 



We will briefly outline the result. The exact solution of the above set of equations 
reads: 



^t{'rV-oj)t 



J{jV-u>)t 



iV tiiii -I Kt\ i sin J Kt ^ , , 
cos jKt a(0) J- <7_(0) 



cos ^Kt + 



K 

iV shi^Kt 

k 



K 

a_(0) ^a(O) 



(22) 



where K = + 1. Note that K and V commute so that their mutual 

ordering is irrelevant. 

We study the case in which the electromagnetic field is in a coherent state, a con- 
dition that coincides with the common experimental situation of a resonant laser 
mode interacting with a spin-i system. Assuming initial factorization between S 
and A, the density matrix of the total system S+A at time t = 0, is: 



1 °° °° p-l"r 

Pit = 0) = -(1 + {a,)a,) Tffr^"'"*"!^)^"^! 

fe=Om=0 VK.VTO. 



(23) 



where |A;) denotes the ket for the photon quanta and lap is the average photon 
number in the coherent state. 

We consider as possible triplets of commuting observables the following: 



t; 



(Tz , a' a , a^a'a 

but others combinations are possible. A direct evaluation shows that only the second 

choice produces a set of independent equations relating the measurements of the 
chosen commutating observables at time f > 0, after turning on the interaction, and 
the state of the system 5 at time t = 0. Only in this case then, a reconstruction of 
the state of the spin-i system at time f = is possible . We report here the details 
of the calculation only for this relevant choice of observables. We have to calculate 
the expectation values of these three observables at time t. In order to perform this 
calculation, we will make use of the eigenfunctions of the V operator, which are: 



l*f) = l« = |0)|-> (" = 0), (24) 



and the operators V and K, when applied to these functions, evaluate to: 

v\ci>t) = ±V^\4>t), k\ci>t) = sf^\ci>t). 



(25) 
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Let us introduce the following notation: 



f{V) 
S{K) 



^i(-fV~uj)t 

i sin ^Kt 

k 



ft = {^i\fiy)\€) = e^^^-'^--)* - {ftr = It 



Sn = {4'n\S{K)\(l}^) = j== 

\/n + l 



g{V, k) = cos-ykt + VS{k) ^ = (0± \g{V , k)\4>^) = cos7\/^t ± ^/^S^ ^ {gtf = 9. 



where, for the sake of compactness, a bar in some cases is used instead of the asterisk 
to indicate the complex conjugation operation. In this notation it holds that: 

a+{t) = {a+g\V,k) - atSt(^))/t(F) (26) 
a^a{t) = {a^g{V, k) - a+S\k)){g\V, k)a - S{k)a.) 
a^aa+{t) = {a''g{V,k) - a+S\k)){g\V ,k)a - S{k)a.){a+g\V ,k) - a^S\k))p{V) 

For a generic observable O we define: 

{d{t)) = Tr[p{t = o)d{t)] = Er=o j:i=±{mt = o)o(i)ic) 

We can then proceed to evaluate the expectation values at a generic time t of the 
chosen triplet of commuting observables. As a first step we have: 



n=0 



fn9n fn 9n / — c ^ 1*^1 ft fn 



2a 



a 



a ft fn ^ -v I f ft 9n ~l~ fn 9n l — a ft fn 

^n ^^2 + t: Vnbn ^ 



Ai- 



- f + - f~ 

\* Q J n J n ^* 
^2 ~ n ^ 



n=0 



n! 



+ \Sn? Ai + 



sJ-^^{aXl-a*X2) + \a\^\g+\^ 



(27) 



n=0 



^ f |^+|2/ rzfn 9n fn 9n q fn 9n fn 9n 9n 9n ■ 

— \\9n\ 7, 7, T 



O ft - fn /„ W?{n + '^) ,, 1 ^^ , o ftgt + fn 9n 9t + 9n W? M 1' 

2 n ~ A^^'' + 2 2 ~ aT' 



+|5„P( 



ft9n - In fin , ft + fn q + ^ 9n + 9n ft " /„ 



?n 9n 



2y/n 2 " y/n 2 2 ' " 2 

ft fn \'^\ 9n ~^ 9n ft ~ fn \\ \ , f ftSn ~^ fn 9n i +|2 , a ft ~ fn 

^—-7^^ — 2— v^^+r — 2 — — 

\g+\^V^ + 2Sn 



ftSn + fn 9n 9n 9n _ |^ |2^2^" ^" ft fn , '1' + 1 



2 _ 2 _ ' _ 2 

9n ~^ 9n f n ~^ f n ■ f n 9n ~^ f n 9n f n fn 



+ 



2 2 

f " + f + 

J n ' J n 



+ ■ 



2n 



+ Sn{\9 



+ |2 



2Vn 



Sn))A^+n-f|5„|2^^±^. 



ft - fn n + 1 f+gt - fn 9n 9n + 9n 



) A2 
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where we have also defined: 



A 



A2 = (a+(0)), A; = (a_(0)) 



(28) 



From Eqs. (|27p it is then easy to derive a system of three equations relating the 
three expectation values: 



(a,(t)) = 25R[(a+(0)], (a^a(i)), (a^aa^it)) = 23? [(a^aa+t))] 



(29) 



evaluated at a generic time t, to the variables {ax{0)), {(jy{0)), (^^(O)) , i.e. to the 
density matrix of the spin-^ system at time i = 0. In order to assess if this system 
of equations has solutions, we have to evaluate the determinant of the matrix M 
made up by the coefficients of {(Tx{0)), ((7^(0)), {az{0)) appearing in the system of 
equations determined by the evaluation of (|29p . In order to achieve this goal, it is 

convenient to define My = X^^o ' Ll"^ — -^^ji("-) ^^'^ then determine the coeffi- 
cients Mij{n). Name Aij(ri) the same type of coefficients but relative to the matrix 
A made up by the coefficients multiplying the variables Ai, A2, A2 in the system of 
equations (|27|) . By replacing the definitions of fn,gi, Sn and then simplifying, one 
gets: 



Aii{n 

Auin 
Ai3{n 
A2i{n 



A22{n 
A-ii{n 



-i iV^ cos 7%/ n + It sin ^^/nt 

a 



a r 



> COS 'y^/nt sin 7\/ n+lt 



) (30) 



a^/rl + T 
e*"* cos7\/n 



sin 7\/n+Tt sin ^\fnt 



It cos 'j^/nt 

(1 + 2n)(l + n - |ap) _(! + „ + jap) cos 2-fy/TrTTt 

2(n+ 1) 
ia cos j\/n + It sin jy/nt 



^23 H = -4^2 (n) 

^33("-) = n-^12("-) 



nAii{n), A32{n) = nAi3{n), 

It is easy to check that these coefficients are related to the coefficients Mijin) in 
the following way : 



M,i(n) = 3?[Aa(n)], Af,2(n) = ^[AM+AUn)], M,s{n) = ^[A,2{n)+A*,{n)] 

(31) 

If we now calculate the determinant A(t) of the matrix M , we obtain: 



00 00 00 _ 

EEE- 

1=0 m=On=0 



3|ar|Q,|2(i+n+m) 

l\m\n\ 



e.^kMu{l)M2j{m)Ahk{n) 



(32) 



We evaluate this determinant numerically. Convergence within the seventh signifi- 
cant figure is reached by keeping in each of the three sums the first 30 terms when 
|ap < 9. For larger values of \a\^ convergence turns out to be much slower. In Fig. 1 
the temporal evolution of the determinant is shown. We see that as |ap increases, 
i.e. for larger average number of photons, the determinant has fluctuations of larger 
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amplitude, so that a time with a large enough determinant can be chosen to solve 
for the initial state of the spin- 5 system. When the matrix elements are determined 
with some experimental uncertainty such a choice is a sensible one that allows a 
more accurate determination of the initial state of the spin-| system and avoids 
cases of ill-conditioned matrix inversion. 




50 t 100 150 200 



Fig. 1. Evolution of the determinant for < t < 200, with the following choice of parameters 
of the model: 7 = 0.1, a; = 0.1. The continuous line refers to the case |ap = 1 the dotted line to 
|a|^ = 4, the dashed line to |a|^ = 9, with a real and positive. 



5. Conclusions 

We have illustrated a procedure which allows to reconstruct the state of a spin-i 
system with a simultaneous measurement of the expectation values of three com- 
muting observables, by coupling the system to an assistant. We have also illustrated 
how the procedure works in the simple case of a spin-i system coupled to a coherent 
laser field. We have shown that in this case, after a proper choice of the commuting 
observables, it is always possible to reconstruct the initial state of the spin-i system. 
A radiation source with a large average number of photons allows to implement the 
procedure in an experimentally reliable condition, i.e. with a large absolute value of 
the determinant of the matrix connecting the expectation values of the commuting 
observables at time t to the initial state of the spin-i system. 
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